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I. INTRODUCTION 



The development of laser and evaporative cooling techniques in the last decades has opened up 
possibilities to cool down alkali atoms to temperatures in the micro- to nanokelvin regime. In this 
regime the control over the external and internal degrees of freedom by external fields allow for 
the experimental realization and in the meantime routine production [? ? ? ] of Bose-Einstein 
condensates . Interacting Bose-Einstein condensates (BECs) exhibit nonlinear excitations such as 
solitons [? ? ], vortices [? ? ] and are the origin of nonlinear matter- wave optics in general. 
Macroscopic as well as microscopic setups [? ? ? ] offer the possibility to coherently manipulate 
matter waves. Moreover Feshbach resonances [? ? ] allow for the control of the interaction prop- 
erties. This enables one to create molecular BECs [? ] and to probe the BEC-BCS crossover [? 
] thereby opening up, together with quantum-phase-transition physics in optical lattices, a strong 
link to condensed-matter physics. 

Rydberg atoms, in contrast to atoms in their electronic ground states, possess extraordinary proper- 
ties. They are of gigantic extension and the highly excited states possess extremely long radiative 
lifetimes. Pictorially speaking the valence electron of a Rydberg atom is in a large loosely bound 
orbit. Consequently Rydberg atoms possess large polarizabilities [? ] and are extremely sensitive 
to external fields. Exposing them to a magnetic field they represent paradigm systems for classical 
and quantum chaos as well as modern semiclassics [? ]. A prominent example are the modulation 
peaks in photoabsorption spectra known as quasi-Landau resonances [? ] being a signature of the 
dominant unstable periodic orbits in the classically chaotic regime. 

There are different experimental techniques to trap long-lived cold Rydberg atoms being of op- 
tic [? ], magnetic [? ] or electric [? ] origin. Specifically for highly excited Rydberg atoms, 
their electronic structure cannot be neglected. Therefore it is crucial to study electronic proper- 
ties of Rydberg atoms in inhomogeneous magnetic-field configurations putting aside the crude 
model of a point- like atom which is adequate for atoms in their electronic ground states. First 
investigations assumed a fixed nucleus taking advantage of the large atomic mass compared to the 
electronic mass. Major results observed were interwoven spin-polarization patterns and magnetic - 
quadrupole-field-induced electric dipole moments induced by magnetic quadrupole fields [? ? ]. 
Taking into account the finite mass of the nucleus (ionic core) leads already in the presence of 
homogeneous magnetic fields to effects due to the non- separability of the center of mass and elec- 
tronic motions [? ? ]. Trapping of Rydberg atoms in the quantum regime, i.e. a quantized center 
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of mass motion coupled to the electronic state manifold, has been demonstrated very recently in 
case of the 3d quadrupole field for high angular momentum state [? ] and for the Ioffe-Pritchard 
configuration [? ]. 

Starting from an ultracold atomic cloud of ground state atoms, ultracold Rydberg gases can be 
created. Theses gases exhibit a new kind of many -body physics [? ? ]. Since the thermal motion 
can be neglected on the corresponding relevant short time scales, one can also speak of frozen 
Rydberg gases. Rydberg-Rydberg interaction, arising from the large polarizability and the large 
dipole moment of the Rydberg atoms, are of major importance for these systems. This interaction 
can be controlled by tuning external fields. In particular the strong dipole-dipole interaction of the 
Rydberg atoms can lead to the so-called dipole blockade which inhibits the creation of a second 
excited atom in the vicinity of a Rydberg atom [???]. Due to these unique properties ultracold 
Rydberg gases represent an intense topical field of research both theoretically [? ] and experimen- 
tally [? ]. 

The problem we investigate in this paper bridges the gap between the situations described above, 
namely the cases of many Rydberg atoms and that of many ground state atoms: We explore the 
spectral structure of a Rydberg atom immersed into a Bose-Einstein condensate, assuming that 
the dimension of the condensate is much larger than the size of the Rydberg atom. This system 
is within the above-outlined context of immediate physical interest and well-accessible within 
present-days ultracold experimental techniques. We start our investigations with the many body 
Hamiltonian. This Hamiltonian is expressed in second quantization with the help of the field op- 
erators. We consequently derive coupled equations of motion for atoms in the electronic ground 
state and in an excited state. Employing several approximations allows us to simplify and decouple 
these equations. Consequently the Gross-Pitaevskii equation appears for the ground state atoms. 
For the excited atom the electronic energy levels depend on the global appearance and the symme- 
try of the density distribution of the ground state atoms. We investigate the cases of homogeneous, 
isotropic harmonic and axial symmetric harmonic density distributions of the ground state atoms 
figuring out a crucial dependence of the electronic spectrum on the shape of the density distribu- 
tions. 

In section II we start with the many-body Hamiltonian of the system and derive the Heisenberg 
equations of motion of the field operators for a Rydberg and a ground state atom, respectively. Sec- 
tion III contains our approximations and their justifications. We introduce contact potentials for the 
interactions and neglect the backaction of the Rydberg atom onto the condensate thus separating 
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the corresponding equations of motion. In Section IV the spectral properties of the Rydberg atom 
for different density distributions of the condensate are investigated. We study a homogeneous 
condensate, an isotropic harmonically and anharmonically confined condensate and a condensate 
with an axially symmetric density distribution. Finally Section V provides the conclusions and a 
brief outlook. 

II. HAMILTONIAN AND EQUATIONS OF MOTION 

In our model we consider atoms consisting of a point-like ion and a valence electron. Let us 
distinguish between ground state atoms with non-excited valence electrons and Rydberg atoms 
with highly excited valence electrons. Assuming that only one atom is excited into a Rydberg 
state we separate the interaction of this highly excited atom with the ground state atoms into two 
parts: The excited electronic and the ionic collisional interaction. This decomposition is good for 
high-angular-momentum Rydberg quantum states. Furthermore we neglect mass correction terms 
due to the finite mass of the core so that the center-of-mass of the atom is located at the position 
of the ion. The many-body Hamiltonian is then given by 

i i e i i i,j 

+ piv IA (R* - Rj)i± + p i veA (Ri, Rj) p i d) 

where M, m e are the atom and the electron mass, respectively. (P i: Rj) denote the center of 
mass momentum and position vector for the i-th atom, respectively, while (p i: r*) represent the 
corresponding momentum and coordinate vector for the valence electron with respect to the %- 
th parent ion. V e denotes the interaction of the ion with the valence electron which accounts 
in principle also for the effect of core scattering of the valence electron, e.g., the quantum defect. 
Later we will restrict ourselves to a pure Coulomb potential which correctly describes high angular 
momentum Rydberg quantum states. The trapping potential V T is based on an adiabatic coupling 
of external fields to the magnetic moment or polarizability of the atom. V AA , V IA and V eA are the 
interaction potentials between two ground state atoms (AA), an ion and a ground state atom (IA) 
and an excited electron and a ground state atom (eA), respectively. The projection operator P° 
ensures that the atom-atom interaction only occurs between two ground state atoms whereas P 1 
ensures that there is only a contribution when a Rydberg atom interacts with a ground state atom. 
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We define x a (r) to be the eigenfunctions of the electronic Hamiltonian H e = + V e (f), i.e. 
the electronic states with the eigenenergies E e a . For highly excited Rydberg atoms with the excited 
electron possessing a large angular momentum, these eigenfunctions are to a good approximation 
given by the hydrogen wave functions. We will focus here exclusively on this case and will label 
the hydrogen wave functions with the quantum numbers (n, /, m) with their usual meanings. The 
Hamiltonian in second quantization is then given by 



a 

a,/3 

+\j dRj dd¥ Q {^{3)V A \R- 3)i>o{R^i>v{d) 
+ J2j dI *J d^],(R)^^)V IA (R-&)^ a (R)^o(E!) 



+ J2 dR dr /^^(^^(^^^^^(^^^X/s^^^)^^) (2) 

a,/3^0 J J J 

where ^> a is the atomic field annihilation operator associated with an atom in the electronic state 

— * — * 

a. Consequently ^ a (R) describes an atom located at the position R in the electronic state \a). In 
the above notation we suppressed the time argument. One obtains via the Heisenberg equation of 
motion for the ground state 

id t MR) = ~^ A MR) + E e MR) + v T (R)MR) 

+ J dR'V AA (R' - R)¥ (R')^ (R')^o(R) 
+ W dR%(R')V IA (R' - R)i>p(R')MR) 
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and for an excited state 

+ E / dr'xl^)V T {R,f) X p{^p{R) 
J 

+ J dR'¥ (R')V IA (R - B!)4f (B!)4f a (R) 

77^0 J J 

III. APPROXIMATIONS : CONTACT INTERACTION AND MEAN FIELD APPROACH 

Since we consider our gas to be ultracold, the velocities of the ground state atoms and of the 
ionic core of the Rydberg atom are very small, such that the potentials V AA and V IA are assumed 
to be determined by a single s-wave scattering length leading to the following parameterizations 

V AA (R-R!) = g5(R-R!) (5) 
V IA {R-R!) = -fS(R-R') (6) 

where g and 7 are the corresponding energy-independent couplings. We assume the interaction 
between the excited electron and a ground state atom to be dominated by s-wave scattering. Being 
interested in principal effects we do neglect at this point higher partial wave scattering which 
might be relevant in specific cases. The velocity of the excited electron cannot be assumed to be 
slow. Therefore we take as a potential for the interaction between the bound excited electron and 
a ground state atom a delta function with a scattering length depending on the kinetic energy of 
the electron. As an example we consider here 87 Rb according to ref. [? ]. A relation between 
the distance of the electron from its parent ion and the velocity of the electron is provided by 
the classical relation E kin (r) = \k 2 (r) = —^2 + \- ^ ne can men determine the coupling as a 
function of the distance r by using the triplet s-wave phase shift 5q given in ref. [? ]. 

r(r) = _ 27r ^fM (7) 

This leads to the potential 

V eA (R,f,R') = T(\r\)6(R + f-Bf) (8) 

For r > 2n 2 where E kin would classically be negative we assume the coupling to be a constant, 
namely the corresponding value for a scattering electron in the zero kinetic energy limit. On the 




E kin [a.u.] 

Figure 1 : The dependency of the coupling F on the kinetic energy of the electron. 

other hand side we assume that the Rydberg electron has an upper bound for its kinetic energy 
such that the available data (see Figured]) on T/5 cover the complete collisional regime. The 
maximal kinetic energy is E^ x = 0.00864 a.u. leading to a minimal distance r min = 116 a.u. 
The probability for an electron in a pure isolated Rydberg state, e.g., n = 30 / = 15 to be located in 
a distance from its parent ion smaller than r min is 0.0002 and can therefore be neglected. Insertion 
of the potentials in eq. © for the time evolution of the ground state leads to 

and in eq. © for the excited state to 

a J 

+ J2f dR'T(\R' - R\)¥ (R')xl(R' - R)x-r(& ~ R)M&)%(R) (10) 

Since we restrict ourselves to the case of a single Rydberg atom or several well isolated Rydberg 
atoms in the not-interacting approximation one can substitute the field operator associated with 
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an excited atom ^/ a (R, t) by a wavefunction ip a (R, t) where a, as already indicated denotes the 
electronic state. It is important to note that the ansatz of a totally symmetric bosonic wave function 
of ground state atoms all residing in the same spatial orbital implies that any of these atoms can be 
excited in the course of the preparation of the system of a Rydberg atom plus condensate. Within 
our approach this property is reflected by the fact that we have an equally delocalized density of 
both the ground state atoms as well as the Rydberg atom (see further below). 
Furthermore we assume the ground state to be occupied macroscopically such that we can intro- 
duce a mean field ip(R, t) for the ground state atoms. This leads to 

id^(R,t) = (^-2_A + E e + V T (R) + g\i;(R,t)\ 2 + 1 J2\M^t)\ 2 

+ J2 I dR!Y{\R-R!\)r p {B!,t) x %R-R!) 
x ,(R-R')^(R',t)~)ij(R,t) (11) 

and for the excited state to 

+ E / dr'xl{7)V T {R^) X p{?)MR,t) 

(3 J 

+ J2 j ^rd^ - ^Dx^C^ 7 - ^)x-y(^ - ^)|V(^ , ,t)lV-y(^,t) (12) 

7^0 J 

Eqs. (II 1112)) represent coupled integro-differential equations for the mean field of the condensate 
and the center of mass amplitudes for the electronic states {a} of the Rydberg atom. Eq. dTTb 
contains the term of the traditional Gross-Pitaevskii equation augmented by terms that couple the 
mean field of the ground state atoms to the Rydberg atom. Eq. (fT2l) is the equation of motion for 
the center of mass amplitudes of the individual states {a} being coupled to the density of the mean 
field. 

Eqs. (|1 1112)) represent a major challenge concerning their numerical solution. The simultaneous 
impact of both the 'background' coherent matter wave of ground state atoms and the trap po- 
tential, originating from electric, magnetic or electromagnetic fields on the Rydberg atom is of 
most delicate nature. Moreover since the focus of the present work is to reveal the impact of the 
matter wave background on the spectral properties of the electronic Rydberg wave function we 
will proceed in the following with a specification of the setup justifying the neglect of the terms 

8 



involving V T in eq (fT2l) . At first a condensate of ground state atoms is created. Then the trapping 
potential necessary for the formation of the condensate is switched off. Since the atoms in the 
condensate are ultracold one can neglect expansion effects of the atomic cloud on the timescale 
of the excitation of an atom into a Rydberg state via a laser. Due to this procedure one needs to 
consider the trapping potential in eq. (fTTj) for the ground state atoms since the trap defines the 
shape of the condensate. However in eq. (TT2l) for the Rydberg atom the trapping potential should 
be absent. Because of the the slow expansion of the condensate cloud the situation after the exci- 
tation process can be regarded as a stationary one. The equation for the stationary solution with 
the chemical potential fi for the ground state atoms then reads 

^(R) = (-^A + E e + V T (R)+g^R) 2 +jJ2MR) 2 

An estimation of the impact of the individual terms on the energy of the condensate with N be- 
ing the number of atoms leads to Eaa ~ gN 2 for the energy resulting from the interaction of the 
ground state atoms with each other, Eja ~ 7^V and E e A ~ TN for the energy contributions result- 
ing from the interactions of the ground state atoms with the ion and the excited electron, respec- 
tively. For a large number of ground state atoms in the condensate one can therefore neglect the 
contribution of V IA and V eA in comparison to the contribution of V AA . The above approximations 
lead to a decoupling of the mean field equation for the condensate from the 'field' of the Rydberg 
atom. At the same time the equation for the Rydberg atom remains coupled to the condensate (see 
below, e.g., eq. (TT61)). Additionally we omit the energy offset due to the electronic energy of the 
ground state atoms. After introducing the normalized wavefunction 4>o(R) = -^=ip{R) and the 
density function p(R) = ip(R) 2 the usual Gross-Pitaevskii (GP) equation appears for the ground 
state atoms [? ]. 

VMR) = (-^& + V T (R)+gp{R)y o (R) (14) 

Since we are in the present work interested in illuminating principal effects of ultracold clouds of 
ground state atoms on Rydberg atoms it is at this point appropriate to firstly employ the Thomas- 
Fermi approximation to the solution of the GP equation. 

P (R) = fi ' V g T(R) Q(R 2 -R 2 ) (15) 
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The stationary equation belonging to eq. (fT2l) reads 

+T(\R'\)p(R' + i?)) X7 (i?> 7 (^) (16) 

where the density p is determined by eqs. (1141151) . In order to proceed we make yet another 
crucial approximation. As discussed above we assume that the condensate wave function re- 
mains unchanged in the course of the switching-off of the trap and the excitation of the Ry- 
dberg atom. Now we additionally assume that the center-of-mass state of the Rydberg atom 
is the same as the one of the mean-field of the condensate. The motivation for this approx- 
imation is the fact that the original totally symmetric microscopic many boson wave function 
of ground state atoms leads to equal probability of Rydberg excitation for all atoms. Conse- 
quently, the excitation process leads to a collective symmetric many -body quantum state \ijj e ) = 
■^W+f l^)il^)2|^)3 • • • \g)i-i\e)i\g) i+ i . . . \g)N+i where \g) t represents an atom numbered % 

in the ground state and |e)j an atom % in a Rydberg state. It is therefore natural to assume that the 
center of mass amplitude for the Rydberg atom equals the corresponding normalized mean-field 
amplitude of the ground state atoms. 

The above concept can be implemented by expanding the center-of-mass components ipa of the 
Rydberg atom in a basis set {(f>i(R)} il} a (R) = J2i 4>i{R)ha with (f) being the macroscopic wave 
function. Then all other expansion coefficients can be neglected compared to the expansion coef- 
ficient of 4> . 

MR) = UR)ha « MR)b 0a (17) 

i 

In the following we will suppress the index denoting the center-of-mass state. Insertion of the 
above approximation and projection on <p = (fi leads to the matrix equation 

J2 H « b -y = ^ (18) 

with 

HI = J dRj d^( J R) X *( J R')(^ + 7P(^)+r(|^ / |)p(^ / + J R))x 7 (^)0(^) (19) 

The first term in the above matrix is the electronic energy. Labelling the hydrogen wavefunctions 
X a with the quantum number (n, I, m) and assuming that V e (f) is a pure Coulomb potential and/or 
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focussing on higher angular momentum states in case of many-electron atoms with core scattering, 
is a constant within a manifold of states belonging to the same principal quantum number n. 
We also focus on the case where the terms providing the mixing of the degenerate hydrogenic 
eigenstates do not induce mixing of states belonging to different n-manifolds. This implies that 
the interaction energy of the Rydberg atom with the ground state atoms is smaller than the spacing 
between these manifolds (see also section IV). The second term results from the interaction of 
the ground state atoms with the ionic core. It does not depend on the Rydberg state and therefore 
cancels if one observes energy differences between two excited states. The last term results from 
the interaction of the excited electron with the ground state atoms and depends strongly on the 
electronic state. Consequently we will focus on this term. The corresponding matrix reads 

Ml = J dRj ^(^)Xa(^)r(| J R'|)p( J R' + J R)x 7 ( J R / )0(^) (20) 

M2 are the matrix elements of the effective potential 

V(R') = J dRct){R)T{\R!\)p{R' + R)<j>{R) (21) 

in hydrogen wavef unctions. This potential can be varied in a wide range by changing the density 
distribution of the condensate due to the external trapping potential. As we shall exemplify in sec- 
tion IV the effective action of the condensate on the Rydberg atom via the potential ((211) can lead, 
for example, to the case of a diamagnetic Rydberg atom if one chooses a strong two dimensional 
confining harmonic trapping potential 

IV. SPECTRA OF RYDBERG ATOMS IN A CONDENSATE 
A. Homogeneous condensate 

We first consider as a trapping potential a rotational symmetric box of radius Rq. 

V T (R) = V (l - &(Rl - ^ 2 )) (22) 
This potential leads, for a large value of Vq, to a density distribution which increases from almost 



zero at the edge of the boxes to a constant on a scale of the healing length £ = y snp a which is 
determined by the scattering length a of the atoms. For a large enough box one can neglect these 
edge regions of non-constant density values and assume the density to be constant within the box 



11 



and zero outside. 

p(R) = Po Q(R 2 - R 2 ) (23) 
This leads with the number of particles being N to the normalized wave-function 4>{R) = y/fj- 



Insertion of the density distribution (|231) in eq. (1201) leads to 

mi = j ' diij d^(fl)x*( J R')poe(^-(fl + fl') a )r(|^|)x 7 ( J R')0(fl) (24) 



The six-dimensional integrals (124)) cannot be evaluated analytically and represent a major numer- 
ical challenge. However, a simple approximation render them tractable. Since the 6-function is a 
two- value piecewise constant function dropping the argument R' in the ©-function just means that 
we assume the Rydberg atom alway to be completely 'covered' by the BEC. This leads to 

mi = j dtij ^(^^(^poe^-^rd^Dx^)^) (25) 

For an arbitrary large box this is obviously true but for a finite trap this is also a good approximation 
since we work in the limit where the typical dimensions of condensates are much larger than the 
size of the Rydberg atoms. After inserting the center-of-mass wavefunction <f)(R) one can carry 
out the R integration by utilizing the normalization relation. We label the hydrogen wave functions 
with the quantum numbers n, I, m. If one neglects n-mixing (see section III) the matrix is already 
diagonal since there are no angular dependent contributions. 

KfiZ = ~J dR'R'Rl^R'mR'l) (26) 

The remaining radial integration was done numerically by a Gaussian quadrature. Figure [2] shows 
the dependence of the shift AE R of a certain state on its angular momentum I for a typical density 
Po = 1.5 • 1CT 11 a.u. This shift is negative for all states since the coupling constant of the electron- 
atom interaction is negative in the region where the spatial probability distribution of the electron 
has a maximum. The absolute value of the energy shift increases with increasing n and decreases 
with increasing /. This can be explained by comparing the spatial probability distribution of the 
excited electron and the dependence of the coupling constant on the distance between the ion 
and the electron. The absolute value of the coupling Y increases with decreasing kinetic energy, 
which corresponds to a larger distance between the electron and the ion. At the same time, the 
expectation value of R increases with decreasing / for fixed n and increases with increasing n 
for fixed /. Therefore the absolute value of the energy shift decreases with increasing I. In the 
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Figure 2: Dependence of the energy shift AE R in atomic units on the angular momentum / for different 
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for po = 1.5 • 10" 11 a.u.=10 
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for the case of a homogeneous condensate. 



latter case one must additionally take into account that the spatial dependence of the coupling 
constant r(— + jj) = — 27r tan( ' 5 ° — is different for different n. The energy shift does not 
depend on m, i.e. states with the same quantum numbers n, I are degenerate. For completeness 
we show in Figure [2] all possible values for the angular momentum, although the strict validity of 
the contact interaction potential in s-wave form is only true for higher angular momenta which are 
also required if one wants to neglect core scattering. 

Let us compare our predicted energy shift with other typical energy scales of a Rydberg atom. The 
natural line widths due to spontaneous emission are typically four orders of magnitude smaller for 
states with a large angular momentum than the typical predicted energy shifts. Line broadening 
due to the interaction of the electron with the ground state atoms can be neglected for ultracold 
atoms since the magnitude of the broadening depends on the induced phase shift which is small for 
ultracold atoms. Relativistic effects such as spin-orbit coupling are about three orders of magnitude 
smaller. The pressure shift a Rydberg atom experiences in a cloud of non-condensed atoms is 
exactly the same as the one obtained above if one considers only s-wave scattering for the collisions 
between the non-condensed atoms and the Rydberg atom. 
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6. Isotropic harmonic and anharmonic confinement 



As a next trapping potential we consider an harmonic isotropic potential 

V T (R) = ^Mu 2 R 2 (27) 

This leads to a density distribution of the ground state atoms according to eq. (031) with Rq = 
and a particle number of 

87r , 2a N a u 

N = TsOl <28) 



Insertion of the density distribution in eq. (|20l) leads to 

Ml = j \rJ dR'<P(R) 2 xURi)r(\R'\) ^ - VT{ ^ + R) )q(R 2 - (R! + R) 2 ) X ,(R')(29) 

Typical trap length scales are much larger than typical extensions of Rydberg atoms and the hy- 
drogen radial wavefunctions vanish quickly for arguments larger than the dimension of the atom 
they describe. Therefore we can simplify the argument of the 6-function Q(Rq — (R' + R) 2 ) — > 
0(i?Q — R 2 ) since either R' -C Rq or R' is larger than the extension of the Rydberg atom so 
that Xa(R') ~ 0. This simplification means again (see subsection IV.A) that the Rydberg atom 
is completely located within the condensate. The center-of-mass function (f)(R) is defined as the 
normalized square root of the density. Insertion of (f)(R) and of the hydrogen wave functions 
Xa(R') = Rn.i(R')Yi,m(Q' \ 4>') makes it possible to fulfill all but one integration in spherical coor- 
dinates analytically. This leads for each n manifold to the diagonal matrix 

<{C = ^J O °dR'R' 2 R 2 n , l (R'nR') 

POD 

/ dR'R' A R 2 nl (R')T(R') (30) 
Jo 



Moo 2 



2g 

This matrix is already diagonal because of the rotationally symmetric potential which leads to a 
conservation of the quantum numbers l,m. Again mixing of states with different quantum numbers 
n need not be taken into account because of their energy gap. 

The energy shift AE R of a state for to — » does, according to eq. (|3Q|) . not coincide with the 
energy shift of the homogeneous condensate in eq. (l26l) . This shows that the energy shift depends 
on the shape of the whole condensate and not only on the local density. The overall shape of a 
harmonically confined condensate is also for small trap frequencies to different from the appear- 
ance of a homogeneous condensate. The dependence of the energy shift AE R on the shape of 
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Figure 3: Dependence of the energy shift AE R on the angular momentum I for different n for \x = 1CT 12 



for (a) u = 1 kHz and (b) u = 10 kHz. 

the condensate results from the derealization of the Rydberg atom within the condensate. Figure 
|3] shows the dependence of AE R on the angular momentum I for different n for [i = 1CT 12 and 
co = 1 kHz or co = 10 kHz. Out of convenience we denote the trapping frequency co in Hz whereas 
lkHz=2.42 • 1CT 17 a.u. For to — 1 kHz the first term in eq. (|30l) dominates the energy shift. This 
summand is, apart from a factor 4/7, equal to the energy shift of the homogeneous condensate. For 
to = 10 kHz the magnitude of the energy shift is, compared to 1 kHz, smaller for sufficiently large 
quantum numbers n. For large quantum numbers n the second summand in eq. (|3Q|) becomes im- 
portant since the expectation value of R 2 increases according to the fourth power of n. In Figure 
|4] the dependence of the absolute value of the energy shift for the complete n = 30 spectrum is 
shown as a function of the chemical potential fi for co = 1 kHz. The absolute value of the energy 
shift increases for all states linearly with the chemical potential. This dependence results from the 
linear dependence of the Thomas-Fermi density distribution ([TBI) on the chemical potential. The 
slopes of the straight lines that connect | AEr\ for states with the same angular momentum / devi- 
ate only little. In Figure |5] the dependence of the energy shift for the complete n = 30 spectrum is 
shown as a function of the trapping frequency co for /i = 10~ n . The absolute value of the energy 
shift decreases for all states quadratically with the trapping frequency. The quadratic decrease can 
be explained by the quadratic decrease of the density with the trapping potential (see eqs. (1291301) ). 
The decrease is different for the different states. States with small angular momenta depend more 
strongly on the trap frequency than states with higher angular momenta. 

For an anharmonic trapping potential V T (R) = aR i a similar argument as for the isotropic har- 
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Figure 4: Absolute value of the energy shift \AE R \ of the Rydberg states with n = 30 for an harmonically 
confined condensate with a trap frequency uj = 1 kHz depending on the chemical potential \x. 

monic trap holds thereby leading to the matrix 



Since the potential is rotationally symmetric the resulting matrix is diagonal, too. For a weak 
confining potential the energy shift is apart from a factor 8/11 approximately the same as for the 
homogeneous case. For a stronger confining potential the second term is small compared to the 
third one since the expectation value of R 6 is for typical quantum numbers n, I of Rydberg atoms 
much larger than the corresponding expectation value of i? 4 . As a result one observes a linear 
decrease of the absolute value of the energy shift with increasing confining parameter a similar to 
the harmonic case. 
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Figure 5: Energy shift of the Rydberg states with n = 30 for an harmonically confined condensate with a 
chemical potential fi = 1CT 11 in dependence on the trap frequency to. 

C. Axially symmetric harmonic confinement 

Finally we investigate an axially symmetric harmonic trapping potential 



V T (R) = l -Mu\x 2 + y 2 ) + -MQ 2 z 2 



The density profile is given by eq. (fT51) with R\{&) 



2jx 

Mu 2 (1+(^|-1) cos 2 (6»)) 



(32) 

whereas 9 is the polar 



angle in spherical coordinates. Spatial integration of the density leads to the number of particles 

N = -^-(^i-^Y- (33) 



15g y Mu 2J 

This leads to the analogue of eq. (|29l) with the corresponding density and the center-of-mass 
wave-function defined as its normalized square root. Here we can make a similar approximation 
to the isotropic confined condensate. We assume that the trap dimension in each direction is much 
larger than the size of the Rydberg atom and therefore the Rydberg atom is fully located within the 
condensate. Consequently the argument of the ©-function can be simplified and the integration of 
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the center-of-mass coordinate can be accomplished analytically leading to 

-I roc pi /*27T 

Ml = - dR'R' 2 / dcos9' / ,9' \(t>')V{R') Xl {R' \& ',(/>') 
9 Jo J-i Jo 

x _ - ^Mcos 2 ^ /2 (fi 2 - cu 2 )) (34) 

This matrix can be understood as an effective potential for the excited electron 

V(R>) = E^(^_^!#2_ Imcos 2 ^' 2 ^ 2 -^ 2 )) (35) 

The following commutational and anticommutational relations hold with P x > being the parity 
operator concerning X' and L z > the angular momentum operator concerning the Z' axis. 

{Px>,L z ,} = (36) 
[V,L Z ,] = (37) 
[V,P X '\ = (38) 

If \E, m) is then an energy eigenstate and at the same time an eigenstate to L z > with L z > \E, m) = 
m\E, m) then we have 

L z ,P x >\E,m) = -m\E,-m) (39) 

Therefore \E,m) and Px>\E,m) = \E, —m) are two degenerate energy eigenstates. So one ex- 
pects for all states but states with m = a twofold degeneracy. Doing the remaining angular 
integration with the help of the recursion relations of the Legendre Polynomials leads to the fol- 
lowing non-zero matrix elements within a n manifold 



<S = \~ n r dR'^Rl^RW) 
1 9 Jo 



Mcu 2 ( (n 2 -uj 2 ) Al-m)(l + m) (I - m + 1)(Z + m + 1) 
^ ~a wo? i ivo; TT 



2# V cj 2 v (2/! + l)(2Z- 1) (2/ + l)(2Z + 3) 

r-oo 

x / dR'R^Rn^TiR^R^R!) (40) 
Jo 



M(tt 2 - uj 2 ) 1(1 + m + 2) (I + m + - m + 1)(/ - m + 2) 
2^ y (2/ + l)(2/ + 3) 2 (2/ + 5) 

/•oo 

x / dR'R' 4 R 2 nl (R')T(R') (41) 
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j\ / rn,l—2,m 
n.l.m 



M(Vl 2 -uj 2 ) (l - m)(l - m - 1)(7 + m)(l + m - 1) 



x 



2g V (2/ + l)(2/-l) 2 (2/-3) 

dR'R^R 2 nl (R')T(R') 



(42) 



For n = to the isotropic case is reproduced but in general the matrix is no longer diagonal. 
Because of the axial symmetry m remains a good quantum number. Only in the case of a near 
spherical symmetric trap the angular momentum / represents an approximate constant of motion. 
In general their occurs a mixing of states with even and odd angular momentum I. As expected 
the energy shift AE R of a state depends on \m\, i.e. states where the absolute value of m is 
different are no longer degenerate. We focus on states with n = 30 and restrict our observations 
to a condensate with fixed chemical potential ji = 1CT 11 since the dependence of the energy 
shift on the chemical potential is similar to the isotropic case. In Figure |6] the dependence of the 
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Figure 6: Dependence of the energy shift AE R for n = 30 states on the transversal trap frequency Q for 
u = 1 kHz and pt = 10~ n 



energy shift on the transversal trap frequency VL for u — 1 kHz is shown. The degeneracy of the 
energy shift is reduced to a twofold one. The energy shifts for states with small angular momenta 
accumulate at small absolute values. Therefore / is not a good quantum number. The energy shifts 
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for high angular momentum states remain separated even for a large difference of the longitudinal 
and the transversal trap frequency. Therefore we will concentrate our investigations on states with 
n = 30 and I = 29. 




Figure 7: Dependence of the energy shift of the n = 30 I = 29 states on (a) the transversal trap frequency 
Q, for oj = 50 kHz and (b) the longitudinal trap frequency u> for $7 = 50 kHz with a chemical potential 

H = io- n 

Figure [7J (a) shows the dependence of the energy shift for the n = 30 I = 29 states on the 
transversal trap frequency. The energy of the m = 29 state is almost independent on the transversal 
trap frequency. This results from the fact that the states with extremal m are located predominantly 
close to the X-F-plane. Therefore the dependence of the energy shift on the transversal trap 
frequency decreases with increasing absolute value of the quantum number m of the state. Figure 
[7J (b) shows the dependence of the energy shift of the n = 30 I = 29 states on the longitudinal 
trap frequency. The absolute value of the energy shift decreases for all states since to defines the 
confinement in two dimensions, quadratically with u>. Here the absolute value of the energy shift 
for the m = state decreases more slowly than the corresponding shift for larger m since m = 
states are located around the Z-axis. 

Finally Figure [8] shows the dependence of the energy shift of the m = states. The m = 
states are for small trapping frequencies uj not degenerate since I is no good quantum number. For 
large trapping frequencies uj one recognizes in Figure [8] a rearrangement of the energy spectrum 
leading to a pairing of the energy levels. Pairwise states with different parity which consist of a 
mixing of either exclusively even or exclusively odd angular momentum states experience almost 
the same energy shift. This can be explained by taking the limit of large trapping frequencies in 
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Figure 8: Dependence of the energy shift of the n = 30 m = states on the longitudinal trap frequency u 
for O = 1 kHz and p = 10 -11 



the effective potential (1351) leading to 

V(R') 



Mu 2 T(R!) 
2^ 



R' 2 sin 2 ff 



(43) 



The coupling coefficient T(R') is constant for large distances R' since the electronic energy is close 
to zero. The support of the square of the radial part of the hydrogen wave function times R' A is for 
n = 30 indeed in the regime where T is constant as can be estimated by the expectation value of R' 2 
which is of order of n 4 . In this limit V = const, and for large cu the effective potential reduces to 
the diamagnetic interaction of an atom in a homogeneous magnetic field [? ? ]. Correspondingly 
the resulting spectral structure observed in Figure [8] possess a completely analogous counterpart. 
The fact that a condensate Rydberg interaction can lead to effective potentials that are equally 
encountered for Rydberg atoms in external fields and/or mesoscopic environments is a intriguing 
analogy. 



V. CONCLUSIONS 



We have examined the electronic spectrum of a Rydberg atom immersed in a Bose-Einstein 
condensate of ground state atoms. Starting with the many body Hamiltonian, we expressed the 
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Hamiltonian in second quantization and derive the Heisenberg equations of motions for the ground 
state atoms and an excited atom. For the mutual interaction of ground state atoms and the inter- 
action of the ionic core of the Rydberg atom and the ground state atoms we introduced contact 
interactions with constant couplings, respectively. The interaction of the Rydberg electron and 
the ground state atoms is modeled by s-wave scattering including the dependence of the coupling 
coefficient on the kinetic energy of the electron. The assumption of the macroscopic occupation 
of one state allowed us to introduce a mean field for the ground state atoms. Furthermore we 
neglected the backaction of the Rydberg atom on the ground state atoms, thereby decoupling the 
equation of motion for the ground state atoms from the equation for the Rydberg atom. Thus 
the Gross-Pitaevskii equation appeared for the ground state atoms. The equation for the Rydberg 
atom depends on the density distribution of the ground state atoms. A Rydberg atom in electric, 
magnetic or electromagnetic fields experiences a complicated trapping potential thereby coupling 
its center of mass and its electronic coordinates. In order to avoid the corresponding complica- 
tions we assumed that the excitation of the Rydberg atom takes place after the trapping potential is 
switched-off. This is possible since the expansion velocity of a Bose-Einstein condensate without 
a confining potential is small so that the condensate can be regarded as static on a typical timescale 
it needs to excite an atom. We additionally assumed that the center-of-mass state of the Rydberg 
atom is the same as the one of the mean-field condensate. This is motivated by the fact that the 
original totally symmetric microscopic many boson wave function of ground state atoms leads to 
equal probability of Rydberg excitation for all atoms which suggests to use the same center of 
mass amplitude for the Rydberg atom as for the ground state atoms. This leads to a matrix which 
is analogue to the Hamiltonian matrix of a Rydberg atom in a potential. In our case the effective 
potential is provided by the density distribution of the condensate. We investigated the energy 
shift of a Rydberg state for different density distribution of the condensate. The energy shift is 
always negative due to the coupling constant of the electron-atom interaction being predominantly 
negative in regions where the spatial probability distribution of the excited electron has a max- 
imum. For isotropic condensates the matrix describing the energy shift is already diagonal for 
each n-manifold since the quantum numbers I and m are conserved. Furthermore the energy shift 
does not depend on the quantum number m. The energy shift depends even for an homogeneous 
condensate on the quantum numbers n and I of the state of the Rydberg atom. Its absolute value 
increases with increasing quantum number n and decreases with increasing quantum number 
For an isotropic harmonically confined potential we observed a linear increase of the absolute 
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value of the energy shift with the chemical potential and a quadratic decrease with the trapping 
frequency of the condensate. For an axially symmetric shaped condensate the matrix is no longer 
diagonal and the degeneracy of the energy shifts is lifted apart from a twofold degeneracy for states 
with m 7^ 0. Furthermore the dependence of the energy shifts on the transversal and longitudi- 
nal trapping frequency is crucially different for different states. For a large transversal frequency 
the effective potential corresponds to the diamagnetic interaction of an atom in an homogeneous 
magnetic field. It is amazing that a condensate Rydberg interaction can lead to effective potentials 
analogous to potentials occurring for Rydberg atoms in external fields or more general to Rydberg 
atoms in mesoscopic environments. 
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